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We investigate the effect of the color-flavor locking pairing pattern on the adiabatic radial oscil- 
lations of pure self-bound quark stars using an equation of state in the framework of the MIT Bag 
model. We integrate the equations of relativistic radial oscillations to determine the fundamental 
and the first excited oscillation modes for several parameterizations of the equation of state. For 
low mass stars we find that the period of the fundamental mode is typically ~ 0.1 ms and has a 
small dependence on the parameters of the equation of state. For large mass stars the effect of 
color- flavor locking is related to the rise of the maximum mass with increasing A. As for unpaired 
quark stars, the period of the fundamental mode becomes divergent at the maximum mass but now 
the divergence is shifted to large masses for large values of the pairing gap A. As a consequence, the 
oscillation period is strongly affected by color superconductivity for stars with M > 1.5 M©. We 
fit the period of the fundamental mode with appropriate analytical functions of the gravitational 
redshift of the star and the pairing gap A. We further discuss the excitation and damping of the 
modes and their potential detectability during violent transient phenomena. 



I. INTRODUCTION 

Neutron stars' oscillations have been investigated for 
over almost 50 years (since the work of Chandrasekhar 
in 1964 [l|) using various equations of state (EOS) for 
nuclear matter and quark matter. The first exhaustive 
compilation of radial modes for various zero temperature 
EOSs was presented by Glass & Lindblom in 1983 0. 
Later studies by Vath & Chanmugan Q complemented 
and corrected earlier studies and also considered the case 
of strange quark stars. Further work by Kokkotas & 
Ruoff [J] presented a new survey of the radial oscillation 
modes of neutron stars. Their study includes an exten- 
sive list of frequencies for the most common equations 
of state. Radial pulsations were also considered in the 
case of newborn neutron stars (protoneutron stars) by 
Gondek et al. They found that the spectrum of the 
lowest modes of radial pulsations of proto-neutron stars 
is quite different from that of cold neutron stars. Gen- 
erally, protoneutron stars are significantly softer with re- 
spect to the radial pulsations, than cold neutron stars, 
and this difference increases for higher modes and for 
lower stellar masses [f| . These differences stem from dif- 
ferent structure of proto-neutron stars, which in contrast 
to cold neutron stars have extended envelopes, inflated 
by thermal and trapped neutrinos effects. 

Radial oscillations were also studied in the case of 
quark stars. Stars containing quark phases fall into two 
main classes: hybrid stars (where quark matter is re- 
stricted to the core) and strange stars (made up com- 
pletely by quark matter). It is expected that both kinds 
of stars cannot exist simultaneously in Nature, but it is 
not know which one would be realized (if any). This de- 
pends on whether the energy per baryon of /3-equilibrated 
quark matter at zero pressure and zero temperature is 
less than the neutron mass (the so called "absolute sta- 
bility" condition 0,0])- Analysis made within the MIT 
bag model shows that there is a room in the parameter 
space for the existence of strange stars. Moreover, color 



superconductivity enlarges substantially the region of the 
parameter space where /3-stable quark matter has an en- 
ergy per baryon smaller than the neutron mass @, • As 
a consequence, a "color superconducting strange matter" 
is allowed for the same parameters that would otherwise 
produce unbound strange matter Q . On the other hand, 
within a Nambu- Jona-Lasinio (NJL) description of quark 
matter, the strange matter hypothesis is not favored, at 
least for the most accepted parameterizations of the EOS 
[13]. Thus, stars containing quark phases are believed to 
be hybrid stars within the NJL model. 

Several works have dealt with the problem of radial 
oscillations of strange and hybrid stars [H, IIiTfl6j . Vath 
& Chanmugan Q showed that the oscillation frequencies 
of strange stars have qualitatively a different dependence 
on the central density compared with the case of neutron 
stars, as the periods of all modes go to zero when the 
central density of the strange star approaches its small- 
est possible value. Benvenuto & Horvath [TTJ] also pre- 
sented calculations of radial oscillations of homogeneous 
strange stars, using a parameterized form of the equation 
of state. They showed that the particular form of the 
equation of state allows some simple and general scaling 
relations which may prove to be very useful for the search 
of these objects. Strange stars have been mostly studied 
in the framework of the MIT bag model, but some work 
has been made using other microphysical description of 
quark matter. In particular, Benvenuto and Lugones 
[12 | studied the radial oscillations of strange stars in the 
Quark-Mass Density-Dependent model [II [3 (QMDD 
model). Their results showed that oscillation periods are 
similar to those obtained within the MIT bag model. Re- 
cently, Anand et al. [l3| studied the radial oscillations of 
magnetized strange stars in the QMDD model showing 
that the squares of the frequencies are always decreasing 
functions of the central density of the strange star. In a 
similar way Singh et al. (lij investigated the radial oscil- 
lations of rotating strange stars in strong magnetic fields 
in the QMDD model. They showed that the difference 
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in frequency between rotating and non-rotating stars is 
larger for higher magnetic fields. The change is small 
for low mass stars but it increases with the mass of the 
star. This change of frequency is significant for the most 
massive stars whereas it is marginal for a 1.4 M Q star. 

In the case of hybrid stars Gupta et al. [7J| studied the 
effect of a mixed quark-nuclear matter core on radial os- 
cillations. They found that the effect of the mixed phase 
is to decrease the maximum mass of the stable neutron 
star and to cause a kink in radial oscillation frequen- 
cies at the onset of the mixed phase. This kink can be 
traced to a slight kink in the density profile as well as in 
the EOS for neutron stars having a mixed quark-hadron 
phase [TJj]. Their results appears to be robust and inde- 
pendent of the EOS. Gondek and Zdunik [7J| also stud- 
ied radial pulsations of neutron stars and strange quark 
stars with a nuclear crust. They used neutron star mod- 
els constructed using a realistic equation of state of dense 
matter and strange star models using a phcnomenological 
Bag model of quark matter. They calculated the eigen- 
frequencies of the three lowest modes of linear, adiabatic 
pulsations and found an avoided crossing phenomenon 
that is strongly related to the changes of compressibility 
of the matter throughout the star. 

In this work we shall focus on the radial oscillations 
of strange quark stars paying particular attention to the 
effect of color superconductivity, which for the best of 
our knowledge, has not been taken into account yet in 
the literature. In particular, we shall consider the effect 
of color-flavor locking within the frame of the MIT bag 
model. The paper is organized as follows: in Sec. II 
we present the equation of state of the CFL phase and 
find simple analytical expressions for the relevant ther- 
modynamical quantities. In Sec. Ill we study the radial 
oscillations of CFL strange stars and we calculate the 
fundamental and first excited modes for different values 
of the parameters of the EOS. In Sec. IV we discuss our 
results. 



II. THERMODYNAMICS OF THE CFL PHASE 

As mentioned before, we work here within the hypothe- 
sis that quark matter is absolutely stable and thus quark 
stars are completely made up by an almost symmetric 
mixture of up, down and strange quarks. At sufficiently 
large densities and low temperatures quark matter is a 
color superconductor, which is a degenerate Fermi gas of 
quarks with a condensate of Cooper pairs near the Fermi 
surface (l9j . Color superconducting quark matter can 
come in a multiplicity of different phases, based on dif- 
ferent pairing patterns of the quarks. At asymptotically 
large densities, where the quark masses are negligibly 
small compared to the quark chemical potential, three- 
flavor quark matter is in the color-flavor locked (CFL) 
state [201 ] . In this state quarks form Cooper pairs of dif- 
ferent color and flavor where all quarks have the same 
Fermi momentum and electrons cannot be present [2lj . 



Color-flavor locking has a profound effect on the proper- 
ties of quark matter, mainly on transport properties such 
as mean free paths, conductivities and viscosities. Con- 
cerning the equations of state, the effects enter as a term 
of order (A///) 2 which is of a few percent for typical val- 
ues of the color superconducting gap (A ~ — 150 MeV) 
and the baryon chemical potential ~ 300 — 400 MeV). 
However, the effect is proportionally very large in the 
low pressure regime that affects the absolute stability of 
quark matter. Thus, pure self-bound quark-matter stars 
[strange stars) may exist for a wider range of parameters 
of the MIT Bag model equation of state [8| . This affects 
considerably the mass-radius relationship of CFL stars, 
allowing for very large maximum masses [H, [22[ . 

The equation of state for CFL quark matter can be 
obtained in the framework of the MIT bag model. To 
order A 2 , the thermodynamical potential ficFL reads Q 
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being f2f roo the thermodynamical potential of a state of 
unpaired u. d and s quarks in which all them have a 
common Fermi momentum v, with v chosen to minimize 
f2f rco . The binding energy of the diquark condensate is in- 
cluded by subtracting a condensation term proportional 
to A 2 /x 2 where the chemical potential /j, = (fi u + Hd + 
/i s )/3 is related to v through v = 2fi — (p 2 + ml/3) 1 / 2 , 
being m s the mass of the strange quark. Confinement is 
introduced through a phcnomenological vacuum energy 
density or bag constant B. 

From the r2cFL given above we can obtain the following 
expressions for the pressure p and the energy density e 
to order to 2 M: 
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In order to have the EOS in the form e = e(p), we can 
invert Eq. ((2|) to find [i as a function of p 



fi = —3a 
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and we can write e = e(p) from Eqs. ^ and ([3]) 

with a and \i given by Eqs. (jH) and ^ respectively 
If we need p = p(e) we may write 

e AB 3a fi 2 
3 3 tx 1 
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with given by 
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For the oscillation equation (see Sec. Ill) we need the 
adiabatic index T = (e + p)p~ 1 dp/de. From Eq. ([7]) we 
have 

dp 1 6ap dfi 

de 3 7r 2 de 
The derivative dp,/de can be calculated from Eq. @: 
d^/de = 7r 2 (9/x)~ 1 [a 2 + 47r 2 (e — _B)/9]~ 1 ^ 2 . Thus, we 
have 

dp 1 2a/ 1 \ , in . 

* = 3 + t(j^0' (10) 

which allows to write the adiabatic index T as function 
of p using Eq. (J5)) or as a function of e using Eq. ((8|). 

Since the values of B, m s and A are not accurately 
known we shall consider them as free parameters in the 
equation of state. We emphasize that all the values of B, 
m s and A employed in this paper fall inside the stability 
windows presented in Fig. 2 of Ref. 0; i.e. we always 
obtain strange stars when integrating the stellar struc- 
ture equations. Additionally, the parameters satisfy the 
stability condition to 2 < 2p,A given in Ref. [23j . 

In our calculations of the next section wc shall also 
employ the MIT equation of state for a gas of unpaired 
massless quarks, i.e. e = 3p + 4B. Notice that unpaired 
matter is the same that CFL matter with A = only if 
all the quark masses are zero. If m s is not zero, the Fermi 
momenta of u, d and s quarks arc different for unpaired 
quark matter, but they are still the same for CFL matter. 



III. RADIAL PULSATIONS OF CFL STARS 

In this section we shall study the radial oscillations of 
strange quark stars employing the equation of state of the 
previous section. We shall consider the unperturbed star 
to be composed of a perfect fluid, whose stress-energy 
tensor takes the form 



(11) 



The generic background space-time of a static spherical 
star is expressed through the line element 
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The Einstein equations in such a spacetime lead to 
the following set of stellar structure equations (Tolman- 
Oppcnhcimcr-Volkoff equations) 



dp em I ^ p 



dr 

dv 
dr 

dm 
dr 



2 dp 
e dr 

47rr 2 e, 



4irpr 3 



2m . 

1 ,13 

r 



(14) 
(15) 



where m is the gravitational mass inside the radius r. 
The metric function v has the boundary condition 



v{r = R) = In 1 - 



2M 
~R 
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where R is the radius of the star and M its mass. 
With this condition the metric function v(r) will match 
smoothly to the Schwarzschild metric outside the star. 

To obtain the equations that govern radial oscillations, 
both fluid and spacetime variables are perturbed in such 
a way that the spherical symmetry of the background 
body is not violated. These perturbations are inserted 
into the Einstein equations and into the energy, momen- 
tum and baryon number conservation equations and only 
the first-order terms are retained. 

Several forms of the oscillation equation have been pre- 
sented in the literature. The original form presented by 
Chandrasekhar [l| constitutes a Sturm-Lioville problem 
whose solution provides the eigenvalues and the eigen- 
functions for the radial perturbations (see eq. (59) of 
Ref. p}). The second order oscillation equation given 
by Chandrasekhar can be split into two first order equa- 
tions. This has been done by Vath and Chanmugan Q 
who derived a set of first order equations for the quanti- 
ties Ar/r and Ap/p (see eqs. (9) and (10) of Ref. (H). 
More recently, Gondek ct al. [f| obtained a similar set of 
equations but for the relative radial displacement Ar/r 
and the Lagrangian perturbation of the pressure Ap (see 
eqs. (11) and (12) of Ref. j5f). All these sets of equa- 
tions are equivalent. In this work, we adopted the equa- 
tions of Gondek et al. [1[ , because this system is particu- 
larly suitable for numerical applications and the bound- 
ary condition at the star's surface can be obtained by 
purely physical arguments. Another important advan- 
tage of this system of oscillation equations stems from 
the fact that they do not involve any derivatives of the 
adiabatic index, T. Adopting G = c = 1 the system of 
equations is 
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dr (p + e) 

where u> is the eigenfrequency and the quantities £ = 
Ar/r and Ap are assumed to have a harmonic time de- 
pendence oc e lut . 

To solve equations (fT7|) and (fT8|) one needs two bound- 
ary conditions. The condition of regularity at r = 
requires that for r — > the coefficient of the 1 /r term in 
Eq. (TTT]) must vanish [1, H El] . Thus, we have 
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FIG. 1. The period of the fundamental oscillation mode (in milliseconds) as a function of the mass M and the central energy density e c 
of the star. In the upper panels the results are shown for CFL strange stars with A = 100 MeV and different values of the bag constant B 
and the strange quark mass m s . B is given in MeV fm — 3 and m a in MeV. Dashed lines correspond to m 3 = and solid lines to m a =150 
MeV. In the lower panels the results are shown for hadronic stars employing the Bethe- Johnson EOS (BJ), for strange stars without color 
superconductivity (MITnp), and for CFL strange stars with A = 50, 100, 150 MeV. The curve MITnp is given only for m s = 0. In the 
left lower panel we fix the value of B to 60 MeV fm -3 and in the right lower panel to 100 MeV fm~ 3 . Comparing the curve MITnp with 
the CFL curves for A = 50,100, 150 MeV we notice that the effect of color superconductivity is large for stars with M > 1.5 Mq. In 
particular, notice that the curve with A = 50 MeV and m s = 150 MeV is almost coincident with the curve MITnp with m s = (see 
left lower panel). This happens because the constant a in Eq. (4) is very small in the former case and it is zero in the latter. Thus, the 
equation of state given in Eq. (6) is very similar in both cases. If we compare the curve for A = 50 MeV and for MITnp both having the 
same m s we see that the effect of color superconductivity is non negligible for A = 50 MeV. The curve for A = 50 MeV and m s = 150 
MeV is not shown in the right lower panel because the parameters of the equation of state fall outside the stability window given in Ref. 

a. 



Notice that the eigenfunctions can be normalized in order 
to have £(0) = 1. The surface of the star is determined 
by the condition that for r — > R, one has p — > 0. This 
implies that the Lagrangian perturbation in the pressure 
at the surface is zero. Therefore the second boundary 
condition is 

(Ap) 

surface 

= 0. (20) 

In order to numerically solve the oscillation equations 
we proceed as follows. First, we integrate the Tolman- 
Oppcnhcimcr-Volkoff equations for each set of the pa- 
rameters of the equation of state (B, m s and A) in order 
to obtain the coefficients of the oscillation equations for 
a given central pressure. Then we solve the oscillation 
equations by means of the shooting method: we start the 
numerical integration of Eqs. (TT7]) and (fT5)) for a trial 



value of w 2 and a given set of initial values of £(r = 0) 
and Ap(r = 0) which satisfy at the center the bound- 
ary condition given above. The equations arc integrated 
outwards trying to match the boundary condition at the 
star's surface. After each integration, the trial value of 
w 2 is corrected in order to improve the matching of the 
surface boundary condition until the desired precision is 
achieved. The discrete values of uj 2 for which Eq. ([20]) 
is satisfied are the eigenfrequencies of the radial pertur- 
bations. Our code was able to reproduce the results of 
Vath & Chanmugam Q and of Kokkotas & Ruoff [|J . 

Our results are shown in Figs. 1—4 for different val- 
ues of the parameters B, m s and A of the equation of 
state (falling inside the stability windows presented in 
Fig. 2 of Ref. In Figs. 1 and 2 we show the period 
of the fundamental and the first excited modes of CFL 
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FIG. 2. Same as Fig. 1 but for the period of the first excited mode. 



strange stars as a function of the stellar mass M and 
the central energy density e c . The effect of color super- 
conductivity has two aspects. By one hand, for a fixed 
mass of the star the period of the fundamental mode is 
smaller as A increases. On the other hand, the maximum 
mass of the star increases significantly for large values of 
A. As a consequence, the oscillation period is largely af- 
fected by color superconductivity for stars with masses 
near the maximum mass. In Figs. 3 and 4 we show the 
period of the fundamental and the first excited modes as 
a function of the gravitational rcdshift z at the surface 
of the compact star. The effect of A is similar to the 
observed in the plots of r versus M (c.f. Figs. [T] and 
[2]). For comparison, the calculations were also performed 
for hadronic stars described by the Bethe- Johnson EOS 
[24l | and for quark stars without color superconductivity. 
For low mass stars there is a large difference between the 
oscillation periods of hadronic stars and quark stars, as 
already known from previous calculations without color 
superconductivity (see lower panels of Figs. [I] and [2]). 
The difference is also large for high mass stars but this 
is due to the difference between the maximum mass in 
different models. A similar behavior is found in the plots 
of t versus z (see Figs. 3 and 4). 

Notice that while other EOS parameters such as B 
or m s can have similar effects than A on the oscillation 
periods, the effect of color superconductivity is stronger 



for A ~ 100 MeV or larger. In fact, m s and A only enter 
in the equation of state through the parameter a defined 
in Eq. (|4|) and they have opposite effects. However, since 
m 2 is multiplied by 1/6 and A 2 by 2/3, the effect of m s is 
proportionally smaller that the effect of A. In the case of 
B, we notice that the period of the fundamental mode of 
stars with M < 1.5 M Q is rather independent of the value 
of B (and also m s ) but changes up to a factor of ~ 2 in the 
here-considered range of A. It is also interesting to note 
that, for equations of state of the form P = (e — 4B)/3, 
it has been shown that there exists a simple scaling law 
relating the oscillation periods corresponding to different 
values of B. The scaling law is r' n = (B/B 1 ) 1 / 2 x r„ 
where the oscillation period t' u corresponds to B' and r„ 
corresponds to B (see [ll[ for more details). Thus, for 
that equation of state, there is shifting of the periods for 
different values of B. In our case, a similar shifting is 
observed in the left upper panels of Figs. 1 and 2 but it 
is not governed by a simple scaling law like the one given 
above because of the presence the condensation term in 
Eq. 0, which is not a constant but depends on the 
energy density e. 



For applications it is convenient to have analytic fit- 
tings of the results. The period of the fundamental mode 
can be expressed as a function of the redshift z of the 
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FIG. 3. The period for the fundamental oscillation mode as a function of the gravitational redshift z at the surface of the CFL strange 
star. In the left panel the results are shown for CFL strange stars with B = 60 MeV fm -3 , m a = 150 MeV and different values of A. 
We also show the curves for hadronic stars employing the Bethe- Johnson EOS (BJ) and for strange stars without color superconductivity 
(MITnp). In the right panel the results are shown for CFL strange stars with A = 100 MeV and different values of B (dashed lines 
correspond to m s = and solid lines to m s = 150 MeV). 
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FIG. 4. Same as Fig. 3 but for the period of the first excited oscillation mode. 



star and the pairing gap A as: 



IV. DISCUSSION 



T~n=0 



cli + a 2 S + C13Z + a^z 2 + a^z6 
1 + M + b 2 z + b 3 6 2 + b 4 z 2 + b 5 zS ' 



(21) 



where T n =o is in milliseconds and S = A / (100 MeV). 
The coefficients are given in Table 1 for different values of 
the bag constant B. This fitting is better than a 5% in the 
range MeV < A < 150 MeV and 0.1 < z < 0.98 x z max . 
The maximum gravitational redshift of the star z max as 
a function of B and A can be obtained through 

z m ax = c + CiA + c 2 B + c 3 A 2 + c A B 2 + c 5 B A, (22) 

with B in MeV fm~ 3 and A in MeV. The coefficients 
are given in Table 2. This fitting is better than a 5% 
for MeV < A < 150 MeV and 60 MeV fnV 3 < B < 
120 MeV fm" 3 . 



In this work we have presented a study of the radial 
pulsational properties of strange quark stars, paying par- 
ticular attention to the effect of color superconductiv- 
ity. We have shown that the effect of color flavor lock- 
ing into the oscillation periods is different for low mass 
stars (with M < 1.5 M Q ) and for large mass stars (with 
M > 1.5 M ). 

For low mass stars the period of the fundamental mode 
T n= o is typically ~ 0.1 ms. For a fixed M , t„ = o is almost 
independent on the value of B and m s but decreases up 
to a factor of ~ 2 as A goes from to ~ 100 MeV. 

For large mass stars the effect of color flavor locking is 
related to the rise of the maximum mass with increasing 
A. As for unpaired quark stars, t„ = o becomes divergent 
at the maximum mass but now the divergence is shifted 
to large masses for large values of A. As a consequence, 
the oscillation period is strongly affected by color super- 
conductivity for M > 1.5 M Q . 
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B = 60 MeV fm" 3 


B = 80 MeV fm" 3 


B = 100 MeV fm" 3 


B = 120 MeV fm" 3 


0,1 


6.04684 xlO" 2 


5.21480 xlO" 2 


4.59254 xlO" 2 


4.12120 xlO" 2 


0.1 


-4.39575 xlO" 2 


-2.20430 xlO" 2 


-1.86426 xlO" 2 


-1.55666 xlO" 2 


az 


5.06630 xlO" 1 


3.81115 xlO" 1 


3.39918 xlO" 1 


3.09221 xlO" 1 


04 


-9.53275 xlO" 1 


-7.82398 xlO" 1 


-6.88127 xlO" 1 


-6.26093 xlO" 1 


as 


2.10216 xlO" 1 


1.20291 xlO" 1 


8.61653 xlO" 2 


6.49367 xlO" 2 


bi 


-2.85725 xlO" 1 


-1.73731 xlO" 1 


-1.82305 xlO" 1 


-1.77994 xlO" 1 


b 2 


-1.25543 


-1.81158 


-1.85871 


-1.92614 


b 3 


7.81002 xlO" 2 


5.28224 xlO" 2 


4.17669 xlO" 2 


3.44168 xlO" 2 


bi 


-1.59488 


-4.63609 xlO" 1 


-3.34299 xlO" 1 


-1.82102 xlO" 1 


bs 


8.21011 xlO" 1 


5.2234 xlO" 1 


5.08038 xlO" 1 


4.77180 xlO" 1 



TABLE I. Coefficients of the polynomial fitting of the fundamental period as a function of the gravitational redshift of the star 
Zmax and the parameter & (being 6 = A / (100 MeV)), see Eq. (J2T|. This fitting was calculated for B = 60, 80, 100, 120 MeV 
fm" 3 . In all cases we assumed m s = 150 MeV. 



CO 


Cl 
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4.60362 xlO" 1 


3.62306 xlO" 4 


-8.89339 xlO" 5 


7.93274 xlO" 6 


1.74139 xlO" 6 


-5.88166 xlO" 6 



TABLE II. Coefficients of the polynomial fitting of the gravitational redshift of the star z max as a function of the parameters 
B and A. As shown in Eq. we have z max — c + a A + c 2 B + c 3 A 2 + aB 2 + c 5 BA, with B in MeV fm" 3 and A in MeV. 
We assumed m s — 150 MeV. 



It is well know that the periods of radial oscillations 
of strange stars behave very differently from those of 
hadronic stars [||, specially for low mass stars. This dif- 
ference is amplified by color superconductivity because 
the oscillation periods tend to be smaller as the pairing 
gap A increases. Also, as for unpaired strange stars, the 
oscillation periods go to zero when the central energy 
density moves towards the smallest possible value e m in 
for which the pressure is zero (see Eq. ©). This can 
be understood by noticing that low-mass quark stars are 
very well described as non-relativistic constant-density 
spheres, which have oscillation periods proportional to 
the adiabatic index T (for the fundamental mode we have 
= 47rGp(4r - 3)/3 [H). Since the pressure tends to 
zero everywhere for a quark star with e ~ e m im the adi- 
abatic index T = (e + p)p~ 1 dp/de tends to infinity, w 2 
diverges, and the oscillation period tends to zero. This 
behavior is clear in the right upper panels of Figs. 1 and 
2. 

In addition to knowing the spectrum of pulsations for 
a given internal composition, it is important to realize 
whether the resulting modes are able to survive for a suf- 
ficiently long time provided they are excited in any as- 
trophysically realistic situation. In principle, radial pul- 
sations could be excited during a variety of catastrophic 
events taking place during the compact star's life. Star- 
quakes and stellar collisions are potential mechanisms if 
enough energy is transferred to vibrational modes. How- 
ever, it has been shown that radial pulsations of quark 
stars are quickly damped due to the enormous bulk vis- 
cosity of hot unpaired quark matter [26j . The volume 
oscillation forces the system out of chemical equilibrium 



with respect to the non-leptonic process u + d f-> u + s 
and the semi-leptonic processes u + e~ O d + v e and 
u + e" f> s + v e . In general, semi-leptonic processes 
are slower than the non-leptonic one, and the most ef- 
fective damping reaction in unpaired quark matter is 
u + d <-> u + s [25|, [27j|. Since u and d quarks are es- 
sentially masslcss while m s ^150 MeV, the forward rate 
of 11 + d H 11 + s cannot keep equal to the reverse rate 
and the system cannot keep in equilibrium during oscil- 
lations. This leads to irreversible processes and damping 
in unpaired quark matter. For a typical stellar oscillation 
time of 10~ 3 s, high amplitude oscillations are damped in 
fractions of a second due to the non-leptonic process [28| . 
As a consequence, the detection of any signal related to 
radial stellar vibrations looks unlikely if quark matter is 
unpaired. 

Still, in the CFL phase the contributions to the bulk 
viscosity from the above processes are exponentially sup- 
pressed [2^ |. The thermodynamic and hydrodynamic 
properties are rather determined by the masslcss su- 
pcrfluid phonons if and thermally excited light pseudo- 
Nambu-Goldstone bosons. The contribution to bulk vis- 
cosity from phonons alone (ip <-> tp + tp) has been calcu- 
lated in [10] (see also [3l[ for a more recent analysis). For 
T > 0.1 MeV and a typical oscillation period r = 1 ms, 
the resulting transport coefficient C, v is several orders of 
magnitude smaller than for unpaired quark matter. 

However, depending on the poorly known value for 
5m = rriK° — [J-k 1 the dominant contribution to the 
bulk viscosity may come from weak equilibrium processes 
involving the neutral kaon K° and the bosons tp, e.g. 
ip + ip and K + p <R- ip 0, [32j . For oscillations 
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with a timescale of milliseconds and temperatures above 
a few MeV, the bulk viscosity C, K can become larger 
than for unpaired quark matter (see [29[ and references 
therein). However, ( K is much less than that of un- 
paired quark matter for low temperatures. For example, 
for T ~ 0.1 MeV and Sm = 0.1 MeV, is more than 
six orders of magnitude smaller than £ for unpaired mat- 
ter. Since the dissipation time scale is t a £ _1 [HI, [H| 
we expect that radial oscillations of CFL stars would last 
for several seconds and are potentially observable after a 
catastrophic event involving the quark star. 

Additionally, we must keep in mind that the spheri- 
cal symmetry is broken if the star is subject to rotation 
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